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Abstract. The composition and physical properties of 
several local clouds, including the Local Interstellar Cloud 
(LIC) in which the Sun is embedded, are derived from ab- 
sorption features in the UV spectrum of the star e CMa. 
We derive temperatures and densities for three compo- 
nents by combining our interpretations of the ionization 
balance of magnesium and the relative population of C II 
in an excited fine-structure level. We find that for the LIC 
n e = 0.12 ± 0.05 cm~ 3 and T = 7000±1200 K. 
We derive the ionization fractions of hydrogen and discuss 
the ionizing processes. In particular the hydrogen and he- 
lium ionizations in the LIC are compatible with photoion- 
ization by the local EUV radiation fields from the hot stars 
and the cloud interface with the hot gas. 
We confirm the detection of high ionization species : Si in 
is detected in all clouds and C iv in two of them, includ- 
ing the LIC, suggesting the presence of ionized interfaces 
around the local clouds. 

Key words: ISM: solar neighbourhood - ISM: clouds - 
ISM: ionization - ISM: structure - stars: e CMa - ultra- 
violet: ISM 



1. Introduction 

The star e CMa (B2 II, V = 1.50, vaini = 44kms" 1 , 
£ = 239.8, b = -11.3) (Hoffleit & Jaschek 1982) at a 
distance of 130 pc (Perryman et al. 1997) is by far the 
brightest extreme ultraviolet source in the sky (Vallerga et 
al. 1993) and thus the main photoionization source in the 
Solar Neighborhood (Vallerga and Welsh 1995, Vallerga 
1998). This is principally due to the extraordinary empti- 
ness of the line of sight to e CMa, as we describe below. 
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* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Institute, which 
is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5-26555. 



A beneficial aspect of the emptiness and hence the sim- 
plicity of the sight line is the opportunity for us to study 
with an unusual level of detail individual diffuse clouds in 
the local interstellar medium. This is the subject of this 
paper. 

In particular, the star e CMa provides an opportunity 
to observe the absorption spectrum of the Local Interstel- 
lar Cloud (LIC) surrounding our solar system. In most 
cases we obtained a good signal-to-noise ratio for the ab- 
sorption features, thereby allowing us to derive the chem- 
ical and physical properties within the clouds, including 
depletion, temperature, electron density and ionization. 

The temperature of the LIC is usually derived from 
the width of absorption lines by profile fitting. Different 
measurements applied to various elements all point to a 
temperature of around 7000 K, but with a relatively large 
error. The most precise measurements are provided by the 
observation of the HI La line, because the low mass of 
hydrogen offers the best discrimination between thermal 
and turbulent line broadening when compared to the re- 
sults from heavier elements. Using this method, Linsky et 
al. (1995) found T=7000±900K in the directions toward 
Capella and Procyon. 

The electron density in the LIC has been determined in 
the lines of sight towards Sirius (Lallement et al. 1994) and 
e CMa (Gry et al 1995) using the ratio V(Mgn)/V(Mgl), 
towards 5 Cas with the ratio jV(Nal)/iV(Call) (Lallement 
and Ferlet 1997), and with the ratio JV(Cn*)/JV(Cll) to- 
wards Capella (Wood and Linsky 1997) and the white 
dwarf REJ1032+532 (Holberg et al. 1999). All methods 
give results that are around 0.1 cm -3 , which is roughly 
the same order of magnitude as the neutral gas density. 

There are strong indications that the LIC and similar 
clouds in the Local Interstellar Medium (LISM) are partly 
ionized. The observed fractional ionization of hydrogen 
can be explained by the EUV radiation from white dwarf 
and other stars, in particular (3 CMa and e CMa (Vallerga 
1998). Previous studies of the nearby line of sight toward 
(3 CMa with the Goddard High Resolution Spectrograph 
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(GHRS) on the Hubble Space Telescope (Dupin and Gry 
1998) and an independent facility, the Interstellar Medium 
Absorption Profile Spectrograph (IMAPS) (Jenkins, Gry 
and Dupin 2000) have shown that the two main clouds 
in that sight-line present a very high ionization fraction 
of hydrogen, which could be explained by photoionization 
due to the combination of e CMa and (3 CMa if the clouds 
are located close enough to the stars. 

An outstanding problem is that these and other stars 
do not produce enough photons with energies above 
24.6 eV to explain the high fractional ionization of helium 
in the LISM, as shown by the low value of n(Hei)/n(Hi) 
(equal to 0.07 instead of the cosmic ratio of 0.1, which in- 
dicates that helium is more ionized than hydrogen). The 
two main proposals to explain this phenomenon are i) the 
LISM is still recombining from a much more highly ionized 
state produced by a supernova-related energetic event in 
the recent past (Reynolds 1986 ; Frisch & Slavin 1996 ; 
Lyu & Bruhweiler 1996) and ii) the ionization of He is 
maintained by the diffuse EUV radiation emitted by con- 
ductive interfaces between the cloud edges and the hot 
gas filling the "Local Bubble" in which they are embed- 
ded (Slavin 1989 ; Slavin & Frisch 1998). 

In this paper we present GHRS spectra of e CMa 
over limited wavelength intervals between 1190 A and 
1550 A, ones that include the lines of N I, O I, C II and 
Cn*, S II, Sin, Sim, and Civ at a wavelength resolution 
R~ 100 000. We consider also a profile of Oi at 1039 A, 
recorded by IMAPS at R ~ 60 000. We show how they shed 
light on the knowledge of characteristics of the nearby dif- 
fuse clouds such as temperatures, electron densities, abun- 
dances, and degree of ionization. 

2. Observations and data reduction 

Most of the observations presented here have been per- 
formed in late 1996 with the Ech-A grating of GHRS. All 
data were taken with the C'25 Small Science Aperture 
(SSA), the procedure FP-SPLIT = 4 and a substepping 
of 4 samples per diode (for details of the instrumentation, 
see Soderblom et al. 1995). For data processing, we used 
the standard STSDAS procedures working in the IRAF 
environment. We assigned wavelengths from the standard 
calibration tables. An error of ±1 resolution element on 
the wavelength assignment is expected to arise from mag- 
netic drifts. 

The signal-to-noise ratio (S/N) for all Ech A data is 
about 200 when the flux is at the level of the stellar con- 
tinuum. However there is a degradation of signal quality 
for features that appear in the bottoms of strong stellar 
lines. For most interstellar lines the S/N ranges between 
100 and 200, but it is 80 for Cn and Cn* and 30 in the 
extreme case of Si in where the stellar line is the deepest. 

Observations in the far-UV lines were carried out by 
IMAPS when it was operated on the ORFEUS-SPAS II 
mission that flew in late 1996 (Hurwitz, et al. 1998). 
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Fig. 1. Example of the continuum normalization for the 
Sill line at 1260. 4A. The stellar line is fitted by a poly- 
nomial (shown here superimposed on the spectrum) over 
the velocity interval [-30kms _1 , 40kms -1 ]. 



IMAPS is an objective-grating echelle spectrograph that 
was designed to record the spectra of bright, early- type 
stars over the wavelengths from ~ 950 A to ~ 1150 A 
with a high spectral resolution. For more details on the 
instrument see Jenkins et al. (1996). 

The IMAPS spectra were extracted from the echelle 
spectral images using special procedures developed by one 
of us (EBJ) and his collaborators on the IMAPS investi- 
gation team. The S/N obtained for the interstellar lines 
observed by IMAPS is on the order of 25 to 30. 

In both sets of spectra, there is an uncertainty related 
to the background correction for scattered light on the 
echelle format. For the GHRS Ech-A data, the error in 
the background correction is smaller than a few percent 
of the continuum and is only a concern for strong lines 
that almost reach the zero intensity level in their deepest 
points. In the case of a species for which several lines with 
different oscillator strengths are available, such as Sin, 
the uncertainty is eliminated by comparing the different 
profiles, which must have relative strengths that follow 
their oscillator strengths. For very strong lines like Hi Lya 
and Cn 1334A, the zero level is determined by the base 
of the saturated line. The situation for the 1 1302A line 



is reported in Section 3.1 



Stellar lines were fitted with one low-order polynomial 
over a velocity range covering the interstellar components 
(§||)j by considering the points bluewards of Component 4 
up to V=-30 kms -1 or -40 kms -1 , redwards of Compo- 
nent 1 up to V=+40 kms' 1 , as well as the velocity range 
which is not affected by any interstellar absorption be- 
tween Components 2 and 3. The spectra are then normal- 
ized to these stellar lines to create the interstellar profiles 
with a level continuum. An example of the stellar line fit- 
ting is illustrated in Figure for the Si 11 1260A line. For 
strong or moderately strong lines the uncertainty added to 
the column density estimates by this normalization pro- 
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cess is small. However, for very faint lines (Mgi, C n*, S n), 
it is a major source of error. This uncertainty was taken 
into account when we listed the column density results in 
Table |[ The most difficult case is that for S II, where all 
three interstellar lines arc very faint and appear on top of a 
steep slope of a stellar line. In this configuration there is a 
lot of freedom for the location of the synthetic stellar lines 
which can artificially enhance the interstellar line or, al- 
ternatively, make it disappear almost completely. We thus 
use several options for the S II absorption profiles defined 
by an envelope of ±2<r around the fitted continuum. 

We derive the column densities by using the line fitting 
software 'Owens' developed by Martin Lemoine. Each in- 
terstellar absorption component is represented by the con- 
volution of a theoretical Voigt profile with the instrumen- 
tal profile. The instrumental profile for the GHRS Echelle 
data was assumed to be a Gaussian with a FWHM of 
0.92 diodes (Soderblom et al. 1995). An iterative proce- 
dure which minimizes the sum of the squared differences 
between model profiles and the data points allows us to 
determine the most likely column densities of the absorb- 
ing elements iV(cni -2 ), the radial velocity of the cloud 
(kms -1 ) and the velocity dispersion (b- value) (kms -1 ) 
of each interstellar absorption component. The software 
also allows us to fit the lines from several elements simul- 
taneously, which leads directly to a coherent solution for 
all species in terms of velocity, temperature and turbu- 
lent velocity. The wavelengths and /-values are listed in 
Table [1} 

3. Line of sight structure and column densities 

Previous high resolution spectra of UV absorption lines 
for the elements Fen, Mgil and Mgi obtained with the 
GHRS Ech B revealed a simple structure for the interstel- 
lar medium in the line of sight toward e CMa (Gry et al. 
1995) : only three components were detected in all absorp- 
tion lines, and their contributions amounted to an equiv- 
alent hydrogen column density of less than 10 18 cm -2 , a 
value that is by far the lowest amount ever observed in 
the Galactic disk for a line of sight as long as 130 pc. The 
low column densities offer the advantage that most of the 
UV absorption lines are unsaturated and can yield reliable 
column densities. Small, extra velocity components are de- 
tected in some of the strongest lines, as shown by Gry et 
al. (1995), however none of these exceeds an equivalent H 
column density of about 4 10 16 cm~ 2 . 

The three principal components had heliocentric ve- 
locities of 17, 10 and —10 kms^ 1 . In the previous papers, 
they were identified as Components 1, 2 and 3, respec- 
tively. Two of the components have been identified by Gry 
et al (1995) as the two clouds detected in the line of sight 
to Sirius by Lallement et al. (1994) and described further 
by Hebrard et al. (1999). As Sirius is located at a distance 
of only 2.7 pc, these two components must be situated 
very close to the Sun. 



Table 1. Observed atomic lines. 





Element 


wavelength (A) 


/"-value 


IMAPS 


Oi 


1039.230 


9.20 10~ 3 


Ech A 


Sin 


1190.203 


2.31 10~ 2 




Sin 


1190.416 


2.93 10" 1 




Sin 


1193.290 


5.85 10 _1 




Ni 


1199.550 


1.30 10" 1 




Ni 


1200.223 


8.62 10" 2 




Ni 


1200.710 


4.30 10" 2 




Sim 


1206.500 


1.6710° 




Di 


1215.339 


4.17 10 -1 




Hi 


1215.670 


4.17 10 -1 




Sn 


1250.584 


5.43 10" 3 




Sn 


1253.811 


1.09 10~ 2 




Sn 


1259.519 


1.66 10" 2 




Sin 


1260.422 


1.18 10° 




Oi 


1302.168 


5.19 10~ 2 




Sin 


1304.370 


9.17 10 -2 




Cn 


1334.532 


1.28 10" 1 




Cn* 


1335.708 


1.15 10" 1 




ol IV 


loyo. loo 


0. 14 1U 




Si iv 


1402.770 


2.55 10" 1 




Civ 


1548.195 


1.90 lO^ 1 




Civ 


1550.770 


9.48 10" 2 


Ech B 


Mgn 


2803.503 


3.06 10" 1 




Mgn 


2796.352 


6.15 10" 1 




Mgi 


2852.964 


1.83 10° 




Fen 


2344.214 


1.14 10 _1 




Fen 


2382.765 


3.20 10 _1 




Fen 


2586.650 


6.91 10~ 2 



1 Wavelengths and /- 
by Morton, updating 



values are from a private communication 
data from Morton (1991). 



Component 1 is recognized as the Local Interstellar 
Cloud (LIC) in which the Sun is embedded, for which the 
motion has been characterized by Lallement and Bertin 
(1992). Component 2, called the "Blue Cloud," is also in 
front of Sirius but is distinct from the LIC. Another com- 
ponent (Component 4) was detected on the blue side of 
Component 3 in the strongest lines at —19 kms -1 , while 
a fifth component was detected only in the line of Si in at 
—65 km s~ 1 , and confirmed in H I because it mimics a very 
strong D I absorption (too strong to be due to D I with a 
reasonable D /H ratio) , for its H I absorption turns out to 
coinc ide with the expected D I feature from Component 1 
(see § |3.3| ). Finally, a "Component 0" had been introduced 
in the red wing of Component 1 to improve the fit of some 
of the profiles, but our reanalysis with more complete data 
indicates that this component may not be real. Its intro- 
duction in the previous Ech B data analysis was proba- 
bly a consequence of using a slightly distorted line spread 
function. We have chosen to omit the component in our 
more refined analysis. 

We have performed the line fitting for Components 1 to 
4 over all elements available in the Echelle A and Echelle B 
data, as well as for the Oi line recorded by IMAPS. All 
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Fig. 2. Ech-B (first six plots) and new Ech-A spectra of e CMa (R ~ 100 000), plotted in the heliocentric velocity 
scale (kms -1 ). The histogram-style tracings represent the observations, and the smooth solid lines indicate the fits to 
the interstellar absorption profiles. Dotted lines indicate the assumed stellar continua and the individual component 
contributions. The snectra have nreviouslv been divided out bv a stellar absorntion nrofile. a,s exnlained in Section 0. 
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Table 2. Column densities (cm 2 ) of the five components detected toward e CMa. 



Comp. 
V 


1 

17 kms" 1 


2 

10 kms -1 


3 

-10 kms" 1 


4 

-19 kms -1 


5 

—65 kms -1 


Fen 


1.35 ± 0.05 10 12 


5.2±0.510 n 


1.99 ± 0.1 10 11 


- 


- 


Mgn 


3.1 ±0.1 10 12 


1.05 ±0.05 10 12 


1.39 ±0.05 10 12 


9.0 ± 4 10 10 


- 


Mgi 


7 ± 2 10 9 


5 ± 2 10 9 


7 ± 3 10 9 


- 


- 


Sin 


4.52 ±0.2 10 12 


1.85 ± 0.1 10 12 


1.38 ±0.05 10 12 


2.9 ± 2 10 10 


- 


Sim 


2.3±0.2 10 12 


2.0 ± 1.1 10 11 


1.0 ±0.1 10 12 


2.2 ±0.510" 


2.1 ±0.310" 


Cn a 


(1.4 - 2.1) 10 14 


(0.7- 1.3) 10 14 


5.21 ± 0.2 10 13 


7.3 ± 2 10 12 


2.3±0.310 12 


Cn* 


1.3±0.2 10 12 


2.0 ± 1.1 10 11 


7.4 ± 1.0 10 11 






Oi 


1A±H 10 14 


1.2±0.310 14 


1.93 ± 0.2 10 12 






Ni 


1.70 ± 0.05 10 13 


9.8±0.510 12 


< 2.2 10 11 






Sn 


8.6 ±2.1 10 12 


4.9 ± 1.5 10 12 


3.5±0.910 12 






Civ 


1.2±0.310 12 




3.7±0.410 12 






Shi 


< 310 12 


< 310 12 


< 3 10 12 






Si iv 


< 210 11 


< 210 11 


< 210 11 






Nv 


< 2.5 10 11 


< 2.5 10 11 


< 2.5 10 11 







a the upper limits for Components 1 and 2 are determined from N(S n) as explained in Section 



lines of Oi, Ni, Mgi, Mgn, Fen, Sin, Sim, Cn, Cn**, 
and S II were fitted simultaneously with a unique linc-of- 
sight velocity structure and a unique column density that 
was consistent with different lines of the same species. As 
the absolute wavelength calibration had a precision of one 
resolution element (about 3 kms -1 ), an individual veloc- 
ity shift is allowed for each data set, whereas the relative 
velocity of the components had to be the same for all 
lines. The resulting velocity shifts for all elements have a 
dispersion of ±0.94 kms -1 , in agreement with the preci- 
sion of the wavelength calibration. The derived velocity for 
Component 1 is 16.15±1.5 kms -1 , in agreement with its 
identification as the Local Interstellar Cloud. The veloc- 
ity shifts of the other components relative to Component 1 
are known with more precision : —6.92 ± 0.03 kms -1 for 
Component 2, — 27.50±0.04 kms -1 for Component 3, and 
—36.4 ± 0.3 kms -1 for Component 4. 

The comparison of the synthetic profiles with the data 
suggested that the velocity shift between Component 3 
and Component 1 could be slightly different for Si ill from 
what we found for the other elements. We thus decided 
to decouple the Sim velocity from the determinations 
for other species, and we found that indeed the velocity 
shift between Component 3 and Component 1 for Sim is 
—28.48 kms -1 , i.e., 1 kms -1 more than for the other el- 
ements. This velocity shift however is only one-third of a 
resolution element and must be compared to the velocity 
dispersion of material in the clouds which is more than a 
factor of three higher. 

When we derived 6-values that had the best fit to 
the data, we imposed the requirement that various ele- 
ments should have common values for the turbulent and 



Table 3. b- values (velocity dispersion in kms -1 ) derived 
from the fits. They are given for the elements carbon and 
iron, which represent the two extremes of atomic mass. 
The 6-values for the other elements are intermediate be- 
tween the two. 



Comp. 


1 


2 


3 


4 


Cn 


3.9 ±0.3 


3.2 ±0.13 


4.0 ±0.1 


5.5 ±0.1 


Fen 


3.0 ±0.15 


2.12 ±0.17 


2.44 ±0.02 


4.4 ±0.4 



thermal doppler contributions within each velocity com- 
ponent. Differences in the outcomes for b across differ- 
ent elements were therefore governed only by variations 
in atomic mass. Table || shows the b- values for two ele- 
ments that represent the two extremes in mass: carbon 
and iron. 

The derived column densities for Components 1 to 4 
are listed in Table ||, and the fits of the lines are presented 
in Figure |J Specific considerations for a few elements are 
discussed individually below. 

3.1. Oi column density 

Guidance on the correct zero level of the O I line at 1302A 
was available from the observed zero level of the nearby 
C II line at 1334A which is known to be heavily saturated. 
The presence of clearly visible component structures at 
the bottom of the line indicates that the profile remains 
above zero or, at most, it hits zero over a range of just a few 
points, in accord with the expected noise level. However an 
uncertainty on the exact zero level location still remains, 
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Fig. 3. IMAPS spectra of e CMa (R ~ 60 000) for the 
line of O I at 1039A and Oi** at 104lA. The histogram- 
style tracings represent the observations and the smooth 
solid lines the fits to the interstellar and telluric absorp- 
tion profiles. Also shown in dash-dotted lines are the as- 
sumed stellar continuum and the individual component 
contributions. In the Oi** line, only the telluric compo- 
nent contributes to the profile while the O I line is a blend 
of the interstellar components 1 and 2 and the telluric line 
(the telluric component is the central one). 

and this is manifested as an uncertainty in the column 
densities for Components 1 and 2. It is thus useful to check 
the results using the weaker line of O I at 1039A, which is 
available in the IMAPS data. 

The interstellar O I feature observed with IMAPS is 
unfortunately contaminated by telluric O I absorption, as 
revealed by the presence of a strong, narrow O I** absorp- 
tion line at 1041 A. However we can estimate the expected 
telluric contribution for O I. In the Earth's upper atmo- 
sphere, the fine-structure levels should be populated in 
accord with their statistical weights. With this in mind, 
we performed a fit of the far-UV Oi and Oi** lines to- 
gether with the other species (leaving out the O I 1302A 
line), but with the addition of a telluric component sat- 
isfying the constraint that AT(Ol) = 5x AT(Ol**). The 
resulting fits are illustrated in Figure ||. They give results 
for the O I column densities of Components 1 and 2 that 
agree with the range of values derived from the O I 1302A 
line alone, giving an assurance that the background as- 
signment for that line was not erroneous. Table [2] lists the 
1 column densities which are consistent with both 1 
lines. 

3.2. Cn column density 

There are substantial uncertainties with the C 11 column 
density determinations for Components 1 and 2 because 
the line is so strong. Over the velocity range covered by 
Component 1, the profile shows nearly zero intensity - a 
situation that is consistent with arbitrarily high column 



densities. Therefore, while some lower limit for N (C 11) 
could be gathered from the profile fitting, we derived up- 
per limits from the upper limits for ./V(Sn) multiplied by 
the cosmic abundance ratio, (N(Cu)/N(Su)) cosm i C = 20 
(Anders and Grevesse, 1989) . In effect, S11 and Cn are 
both the dominant ionization states in diffuse neutral gas, 
and even if these two elements are in a partially ion- 
ized medium, they should have about the same fraction 
of atoms elevated to higher stages of ionization (Jenkins, 
Gry and Dupin 2000). While this may be true, we also 
know that carbon is usually more depleted than sulfur in 
the ISM, implying that N(Cu)/N(Sn) is probably lower 
than the cosmic abundance ratio. The upper limits for 
N(C 11) listed in Table || for Components 1 and 2 reflect 
these considerations. For Component 3, the line is not as 
saturated, and we derived both the lower and the upper 
limits from the line fitting, and we note that the upper 
limit we measure is very close to the value we would esti- 
mate from S 11. 
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Fig. 4. The four interstellar components of H 1 and D 1 
fitted on the stellar Lyman a profile of e CMa, represented 
by a fourth order polynomial, the coefficients of which are 
also free parameters in the fit. This presentation illustrates 
the complex superposition of profiles. It is impossible to 
find a unique solution for the H 1 and D 1 column densities, 
partly because of the impossibility to define the stellar 
continuum independently. 



3.3. Hi and Dl Lyman a profiles 

In principle, the low column density of neutral matter 
present in the line of sight toward e CMa is favorable for 
studying D/H because the bottom of the Hi line is not 
wide enough to bury the deuterium line. Indeed, as seen 
in Figure [|, the Dl feature is clearly seen in the wing of 
the H 1 profile. Unfortunately, as already noted by Gry et 
al (1995) with the G160M data, it is very difficult to de- 
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rive H I and D I column densities from the Lyman a line 
because 

i) the H I stellar line is narrow compared to the interstellar 
profile and is thus very difficult to define with any accu- 
racy, 

ii) the D I feature is probably blended with an H I absorp- 
tion from the small component at —65 kms -1 . Although 
this component is very weak and probably mostly ionized, 
for it is detected only in Si in and C n, a neutral hydrogen 
column density of only 10 cm -2 could dominate over a 
deuterium feature if D/H ~ 2 1CP 5 . 

iii) since the mass of deuterium is low, its b- value is high 
and the absorption features from the different components 
are not resolved very well, unlike the cases for the other 
elements. This makes it impossible to distinguish the con- 
tribution of each component to the D I profile, which might 
otherwise allow us to separate a blend of three or four D I 
components from the one Hi high velocity component. 

We plot in Figure || some representative fits to the Ly- 
man a profile simply to illustrate the complexity arising 
from the overlapping individual absorbers. In this exam- 
ple the Lyman a profile has been fitted using informa- 
tion from the other elements, without any additional con- 
straints. This solution is not acceptable because it implies 
an Hi column density for Component 2 which is incom- 
patible with the limit set by Oi or Sn. We have checked 
nevertheless that the H I column densities derived from O I 
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in Section 4.1 produce a synthetic Lyman a profile which 
is compatible with the observed profile. 

3-4- Civ profiles 

Some absorption features can be seen on the EchA spectra 
at the bottom of the two C IV stellar lines (Figure ||) . They 
were already apparent at lower resolution in the G160M 
data (Gry et al. 1995). If we fit a stellar continuum outside 
the range where the contributions from Components 1 and 
3 are expected , the normalized spectra show two compo- 
nents with a velocity separation which is compatible with 
the separation between Components 1 and 3 with column 
densities of 1.2±0.310 12 and 3.7±0.410 12 cm" 2 , respec- 
tively. 

There is a shift of about —5 kms -1 in the absolute 
velocity of the components, compared to the velocities of 
these components measured for the other species. This 
shift is a bit too high to be attributed to an uncertainty 
in the wavelength calibration, which should be less than 

3 kms -1 (this has been verified with the other elements). 3.5. Other high ionization species 



Fig. 5. The Civ profiles of e CMa at 1548.2 A (left) 
and 1550.8 A (right). This fi gure shows the fit of the 
stellar line by a fourth order polynomial (top), followed 
by the residuals arising from interstellar absorption with 
two components (bottom). A fit (solid line) indicates the 
existence of two components with a velocity separation 
equal to that of Components 1 and 3 derived from the 
other elements. However the absolute velocity of ClV is 
shifted by —5 km s -1 compared to the velocities derived 
from the other species. The temperatures implied by the 
fits are equivalent to a kinetic temperature within the 
range 100 000 to 200 000 K. The column densities are given 
in Table | 

inates the possible bias that could arise from the a pri- 
ori choice of the location of the interstellar features when 
defining the stellar continuum. The result of the fit was 
virtually the same as that found with our standard way 
of deriving a normalized spectrum. Once again, a velocity 
separation equivalent to that between Components 1 and 
3 was found. This separate exercise reinforces our proposal 
that the C IV lines are real and not simply a consequence 
of our incorporating information from other elements in 
lower stages of ionization. 



However, velocity discrepancies between highly ionized 
species and lower ionization stages are not unusual; see, 
e.g., Sembach, Savage and Jenkins (1994). 

We have also performed an evaluation of the C IV lines 
with a fit to the stellar and interstellar lines simultane- 
ously over a velocity interval that includes the region cov- 
ered by Components 1 to 3 and without specifying the ve- 
locity of the interstellar components. This approach elim- 



Absorption features from the O vi doublet have been 
looked for in the IMAPS spectrum but they are not vis- 
ible. None of the other high ionization species Siiv, S in 
and N v have been detected in the GHRS spectrum. 



4. Neutral and total gas column densities 
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4--1- Neutral gas column densities 

The ionization fractions of oxygen and nitrogen are 
coupled to that of hydrogen via resonant charge ex- 
change reactions. In particular, the rate coefficient for 
charge exchange of On with Hi is exceptionally strong, 
- lO^cnrV 1 (Field & Steigman 1971), and Oi is there- 
fore a very good tracer of Hi. We can thus derive the 
column densities of neutral hydrogen in the components 
from their respective O I column densities through the re- 
lation JVoi(Hi) = #(Ol) / (0/H)i SM , where (0/H) JSM is 
the abundance of oxygen relative to hydrogen in the in- 
terstellar medium. If we adopt the value derived by Meyer 
et al. (1998), (0/H) ISM =3.16 1CT 4 , wc find that 
iVoi(Hl)=4.4+Jj 10 17 cm- 2 for Component 1, 
AT OI (Hl)=3.8±0.9 10 17 cm- 2 for Component 2 and 
i\T O i(Hl)=6.1±0.6 10 15 cm- 2 for Component 3. 
With these numbers, we arrive at a total Hi column den- 
sity in the range 7 10 17 to 1.1 10 18 cm~ 2 , which is compat- 
ible with that derived from the absorption of the extreme 
ultraviolet flux from the star based on EUVE spectra (Val- 
lerga et al. 1993 ; Cassinelli et al. 1995). 

If we adopt a different tactic by assuming that the 
diffuse clouds in the local ISM are undepleted and the 
oxygen abundance is equal to the abundance in B stars, 
(0/H)*=4.68 10~ 4 , then the Hi column density is lower, 
ranging from 4.8 10 17 to 7.2 10 17 cm -2 for the whole line 
of sight, with 

7V OI (Hi)=3.0±o 4 10 17 cm~ 2 for Component 1, 
AT OI (Hi)=2.6±0.6 10 17 cm~ 2 for Component 2 and 
iV O i(Hi)=4.1±0.4 10 15 cm- 2 for Component 3. 

We can also derive the H I column densities implied by 
the Nl column densities. If we again use the abundance of 
nitrogen in the interstellar medium derived by Meyer et 
al. (1997), (N/H) /S A/=7.51(r 5 , we find 
iV Nl (Hi)=2.1±0.210 17 cm~ 2 for Component 1, 
7V Nl (Hi)=1.3±0.1 10 17 cm~ 2 for Component 2, and 
7V Nl (Hi)<2.910 15 cm~ 2 for Component 3. 

We note that in all three components, the neutral col- 
umn density derived from N I is significantly lower than 
that derived from O I. This N I deficiency in the local ISM 
has already been noticed by Jenkins et al (2000) from 
far-UV spectra of white dwarf stars observed with FUSE. 
They have shown that the Nl deficiency, as that of An, 
favors the existence of a source of ionizing photons with 
E>24.6 eV in the local ISM to explain the Hei ionization. 
For the subsequent discussions, we will adopt for N(Hi) 
the ranges derived from N(Oi), keeping both options for 
the oxygen depiction. 



4-2. Total gas column densities 

S II is traditionally used as an indicator of the total (i.e. 
neutral plus ionized) gas column density in the interstel- 
lar medium because i) sulphur has little or no depletion 
onto dust grains (see, e.g., Savage & Sembach (1996) and 



Fitzpatrick & Spitzer (1997)) and ii) the ionization po- 
tential of S II (23 eV) is high and therefore S II is of- 
ten assumed to be the dominant ionization stage both 
in HI and HII regions. If we trust that assumptions i) 
and ii) are valid, we can derive the total column den- 
sities for Components 1, 2 and 3 from the expression : 
Ng n (H tot ) = N(Sn) / (S/B.) cosmic, where (S/H) cosmic is 
the cosmic abundance ratio taken from Anders & Grevesse 
(1989) : (S/H) cosmic =1.8410- 5 . In this context, we obtain 
N Sn (H tot )=4.7± 1.1 10 17 cm- 2 for Component 1, 
AT SlI (H tot )=2.7±0.8 10 17 cm- 2 for Component 2, and 
A Sn (H tot )=2.0±0.5 10 17 cm- 2 for Component 3. 

5. Electron density and temperature 

We follow the method of Jenkins, Gry & Dupin (2000) to 
derive both the temperature and the electron density in a 
component by combining the information from two ratios 
which depend in different manners on these two physical 
parameters. 

The relative populations of the fine-structure levels of 
C II are governed by the balance between collisions and 
radiative de-excitation. In the diffuse warm medium, the 
collisions are dominated by electrons. The condition for 
equilibrium 



n(e)7i )2 n(Cli) = [n(e)72,i + A 2 ,i]n(C II*) 
leads to an equation for the electron density 



n(e) 



92A 2>1 T - 5 


"n(C II*) 
. "(ClI) 




8.63 io- 6 n lia {(^) 


exp (^a) 




"n(C II*)" 
_ n(ClI) _ 


} 



(1) 



(2) 



where the radiative decay probability for the upper level 
is A 2< i = 2.29 lO^s" 1 (Nussbaumer & Storey 1981), and 
the collision strength £1 12 = 2.81 (Hayes & Nussbaumer 
1984). Equation || thus gives a curve of n(e) versus T for 
any value of the ratio n(Cll*)/n(Cn). 

To see the conditions from a different perspective, we 
make use of the equation for the equilibrium between the 
two lowest ionization levels of magnesium which is given 

by 



[r(Mgi) + C(Mgii)n(ff + )]n(Mgi) = 
a(Mgi)n(e)ra(Mgii) . 



(3) 



For the charge exchange rate C (Mg n) that applies to the 
reaction Mgi + H + — » Mgn + H, we used the analytical 
approximation C(Mgn) = 1.74 10~ 9 exp(-2.21 10 4 /T) 
derived by Allan, et al. (1988). We calculate T(Mgi) = 
6.1 10 _11 s _1 at the position of the Sun (for details, see 
Jenkins, Gry & Dupin (2000)). For a(Mgi) we used the 
radiative and dielectronic recombination rates given by 
Shull & van Steenberg (1982), supplemented by the addi- 
tional contributions from low-lying resonance states com- 
puted by Nussbaumer & Storey (1986). As with Equa- 
tion |^, an application of Equation |3| also provides a curve 
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Temperature Temperature Temperature 

Fig. 6. Curves of the electron density n(e) versus temperature T governed by equations |^ and [| for the permitted 
values of the ratios n(Cll*)/n(Cn) and N(Mgn)/N(Mgi). For each of the three components the shaded area where 
the two sets of curves intersect determines the possible values for n(e) and T. 



Table 4. Derived properties of the three main compo- 
nents. 



Comp. 


1 (LIC) 


2 


3 


T (K) 


5700-8200 


8200-30 000 


6000-8400 


n e (cm -3 ) 


0.08-0.17 


0.016-0.088 


0.18-0.28 


n e /n tot 


<0.55 


<0.43 


0.955-0.985 


n tot (cm -3 ) 


>0.14 


>0.034 


0.18-0.29 


p/k (cm- 3 K) 


>1300 


>440 


2300-5000 


length (pc) 


<1.3 


<3.3 


0.17-0.45 



giving n(e) versus T, but this time for any value of the 
ratio 7V(Mgn)/7V(Mgi). 

In Figure ^, we plot for each component the curve n(e) 
versus T of Equation for the maximum and the mini- 
mum values of the ratios n(Cii*)/n(Cn) derived from the 
observations, together with the two curves n(e) versus T 
given by Equation [| for the two extremes for the ratio 
AT(Mgn)/iV(Mgi). The shaded area bounded by the two 
sets of curves defines the possible values of n(e) and T that 
are consistent with both methods for each component. 

We derive 5700 < T < 8200 K with 0.08 < n e < 
0.17cm" 3 for the LIC (Component 1), 8200 < T < 
30 000 K with 0.016 < n e < 0.088 cm" 3 for Component 2, 
and 6000 < T < 8400 K with 0.18 < n e < 0.28 cm- 3 for 
Component 3. 

Note that our results confirm the temperature of 
7000 K commonly found for the LIC in most studies, and 
our electron density of n e = 0.12 ± 0.05 cm -3 confirms 
the ranges found independently in the LIC in other sight- 
lines by Wood & Linsky (1997) and Holberg et al. (1999). 
From our results, the LIC and Component 3 have similar 
temperatures while Component 2 could be warmer. 



6. Ionization fractions and ionization processes 

From the comparison of the total and the neutral hydrogen 
column densities derived from N(S n) and N(0 1) respec- 



tively (see Section ||) , we estimate the ionization fractions 
n e /n tot , listed in Table ^. 

Component 3 is ionized by more than 95 %, making its 
ionization comparable to that of the two main components 
(C and D) in the line of sight to f3 CMa (Jenkins, Gry & 
Dupin 2000). As for Components C and D, the main ion- 
izing source for Component 3 is probably the star e CMa 
and its high hydrogen ionization fraction requires that it 
is located closer to the star, further away from the Sun. 

Component 1 (the LIC) can be at most 55% ionized, its 
maximum ionization fraction corresponding to the case of 
no oxygen depletion in the local cloud. However, formally 
the column density ranges allow the LIC to be neutral. 

Since for Component 2 the neutral gas and total gas 
column density ranges are nearly coincident, Component 2 
is most probably neutral although formally its ionization 
fraction can be as high as 43% in the case of no oxygen 
depletion. 

In order to compare these numbers to the expected ion- 
ization fractions, we have calculated the photoinionization 
equilibria in the local ISM for varying depths of shielding 
by neutral H and He following the equations described 
by Sofia and Jenkins (1998). We consider the effect of 
Vallerga's (1998) composite stellar radiation field supple- 
mented by Slavin's (1989) calculated flux from the cloud 
conductive interface. The results are shown in Figure 0. 
The calculation is the same as that used to produce Fig- 
ure 2 of Jenkins et al. (2000) except for a higher pressure 
p/k, in better agreement with the electron density derived 
in Section || for the LIC as well as with the neutral den- 
sity derived from EUV stellar spectra and inside the helio- 
sphere (Vallerga 1996, Quemarais et al. 1994). We derive 
the expected ionization fraction in the two local clouds 
from Figure ^ and from the estimate of the mean attenua- 
tion of the ionizing flux received in the clouds. Since most 
of the EUV radiation comes from the direction of e CMa 
where nearly all the neutral gas is included in one of these 
two components, the attenuation is derived from their Hi 



column densities given in Section 4.1. 

If we assume that Component 2 is located outside of 
the LIC, further out toward e CMa, it should be shielded 
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Fig. 7. Predicted densities versus H I attenuation column 
density in a local cloud exposed to the local photoioniza- 
tion field (due to hot stars and hot-gas conductive inter- 
face). 

only by its own H I material, with a mean H I attenua- 
tion of 1/2 N(Hi), in all cases lower than 2.3 10 17 cm~ 2 . 
With this upper limit, Figure predicts for Component 2 
a minimum ionization fraction of 0.44, which is inconsis- 
tent with the slightly lower maximum ionization fraction 
permitted by the observations. Note that the temperature 
of Component 2 derived in Section || is higher than that 
used in Figure |?], optimized to match the characteristics 
of the LIC. As the ionization fraction increases with tem- 
perature because the recombination rate is a decreasing 
function of T, the predicted ionization fraction for Compo- 
nent 2 is even higher than that which appears in Figure fj], 
amplifying further the discrepancy with the value derived 
from the observations. This inconsistency would be elim- 
inated if Component 2 were shielded by the LIC, which 
would imply that it is located within the Local Cloud, a 
configuration already proposed by Gry (1996) for some of 
the small components detected close to the Sun in several 
lines of sight. 

For the LIC, the mean attenuation due to the cloud 
itself ranges from 1.5 10 17 to 2.2 10 17 cm -2 depending 
on the assumed oxygen depletion value. From Figure fj], 
this implies an ionization fraction between 0.45 and 0.5. 
If we were to add the attenuation effect due to Compo- 
nent 2 (which is maximum if Component 2 is external and 
shielding the LIC as a whole), the range of possible ab- 
sorbing column is extended up to 4 or 6 10 17 cm -2 and 



we would get an ionization fraction down to 0.34 or 0.24. 
These values are all compatible with the range derived for 
the LIC from the measured column densities. 

Figure [?] also predicts that the ratio n(Hei)/n(Hi) is 
almost constant over the range of Hi column densities 
measured toward the white dwarf stars observed by EUVE 
(around 1 10 18 cm' 2 ) and is close to the measured ratio 
of 0.07 (Dupuis et al. 1995). We conclude from the ad- 
equate prediction of the hydrogen ionization fraction as 
well as of the n(He i)/n(H i) ratio that the photoioniza- 
tion by the EUV radiation field due to the combination of 
the hot stars and the cloud conductive interface is a likely 
representation of the ionization processes in the LIC. 

7. Other physical properties of the clouds 

From the comparison of the ionization fractions and the 
electron density, we infer the total hydrogen density n tot 
(or its lower limit) for each of the three components. These 
are listed in Table ||. In fact the numbers are compatible 
with all three components having a similar total density 
close to 0.2-0.3 cm -3 but different electron densities de- 
pending on their ionization state. 

The thermal pressures in the components are also 
listed in Table |^. To estimate the thermal pressure in the 
clouds from p/k = nT we sum up over all particles : hydro- 
gen (n tot -=nff j+n e ), helium (0.1 n tot ) and electrons (n e ). 
For the LIC, if we adopt the nyj range of Vallerga (1996), 
i.e. 0.15 to 0.34 cm -3 , with our ranges for n e and T, we 
derive p/k = 1900 to 6000 cm -3 K, in agreement with the 
thermal pressure found in Component 3, as is expected if 
they are in equilibrium with a surrounding medium that 
is common to both of them. 

The thicknesses of the components along the line of 
sight are estimated from the ratio of the total column 
density N(H tot ) derived in Section |4~2] to the total density 
n tot . The lengths of the components (listed in Table |J) are 
all very small compared to the length of the line of sight 
: they occupy a total of less than 5 pc. This implies that 
at least 96 % of the sight-line is empty or filled with more 
highly ionized gas. 

8. High ionization species and cloud interfaces 

Substantial amounts of Si in are detected in all compo- 
nents. We find that iV(Sim) = 0.5iV(Sill) in Compo- 
nent 1 and more than 0.7JV(Sill) in Component 3. While 
iV(Si in) is only about 10% the value of iV(Sill) in Com- 
ponent 2, it is 7.5 times AT(Sill) in Component 4. 

In principle Si in cannot come from the same region 
as the other species because charge exchange with even 
small amounts of neutral hydrogen tends to shift the Si to 
lower stages of ionization (Jenkins, Gry & Dupin 2000). 
Nevertheless, it seems clear that Si in arises at velocities 
that coincide with all components, including even the less 
ionized Component 2. This makes it is very likely that Si in 
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is located in regions associated with the clouds, perhaps 
in their outermost layers. 

It is interesting to note that Si III is not detected in the 
spectrum of a CMa (Hebrard et al 1999) up to a limit of 
2 10 11 cm -2 . This limit corresponds to our detection for 
Component 2, but it is more than 10 times lower than the 
column density we derive for the LIC. This discrepancy 
could be interpreted as Si in coming from a completely 
unrelated cloud, which coincidently has about the same 
velocity as the LIC, as proposed by Hebrard et al (1999). 

A difficulty with the above proposal is that this com- 
ponent probably would contaminate other lines, and thus 
it should influence the column density determinations for 
other elements. This hypothetical extra component would 
be likely to have elemental abundance ratios different 
from the LIC -due in particular to different ionization 
fractions- and thus create large perturbations in the col- 
umn density ratios. Yet, if we compare the results for the 
four species for which we have reliable results for both 
a CMa and e CMa sight-lines (i.e. derived from unsatu- 
rated lines) : Ni, Sin, Fell and Mgn, the column density 
ratios between e CMa and a CMa are very similar for 
all elements : N(e CMa)/iV(a CMa) = 1.5 ±0.2c m - 2 . 
As an example of the kind of differences one can expect 
between two different components, the column density ra- 
tios between Component 1 and Component 3 for the three 
ionized species Sin, Fen and Mgn present a dispersion of 
60%. Even worse, the Ni column density ratio between 
Component 1 and Component 3 is more than a factor of 
10 higher than the ratios for the other elements. 

Thus, in view of the small dispersion in the column 
density ratios between the LIC in the e CMa sight-line and 
the LIC in the a CMa sight-line in the various elements, we 
strongly support the idea that there is only one absorbing 
component at the LIC velocity in the e CMa sight-line, 
and thus that the Si III absorption is related to the LIC 
component. 

The fact that N(Si in) is at least a factor of 10 lower in 
the Sirius sight-line than in the e CMa sight-line suggests 
that the region where the Si in originates is extended and 
that more than 90% of it lies beyond Sirius. The possible 
very small velocity shift between the Si in absorption and 
the absorption from less ionized species could be explained 
by the presence of a slight positive velocity gradient to- 
ward the outer layer of the cloud. 

The detection of Sim in the LIC has an interest- 
ing consequence for the relative location of the two local 
clouds. Since by definition the LIC is the cloud in which 
the Sun is embedded, the presence of an extended LIC 
Si in layer past Sirius implies that the LIC and its ex- 
tended layer occupy the full line of sight toward Sirius. It 
follows that Component 2 (the "Blue Component" in the 
a CMa sight-line) should be embedded in the LIC or at 
least in its extended Si in layer, corroborating the sugges- 
tion made in Section pi 



There is a significant difference between the gas re- 
sponsible for the Si in absorption and the gas responsi- 
ble for the C IV absorption : while the widths of the Si in 
profiles are compatible with Si in being at the same tem- 
perature as the less ionized elements, the profiles for C iv 
are clearly broader, implying temperatures of the order of 
100 000 to 200 000 K. This favours the existence of colli- 
sional ionization due to a high temperature. Indeed, our 
C IV column densities and our Si iv upper limits are com- 
patible with the outcome from Slavin's (1989) model for 
the conduction layer between the Local Cloud and the hot 
gas that is supposed to fill the Local Bubble. He calculated 
that AT(Civ) = 2.7 10 12 c m - 2 and iV(Sirv) = 1 10 n c m - 2 
through an interface of this sort, assuming minimal in- 
hibiting effects from magnetic fields. His prediction for 
iV(Siin) of 5 10 10 cm -2 is far below our observed column 
density. 

9. Summary 

We have analysed the interstellar absorption lines in the 
high spectral resolution (R ~ 100 000) UV spectrum of 
e CMa (130 pc). We derive column densities for 11 dif- 
ferent elements in the three main clouds and for a few 
elements with the strongest lines in two additional very 
weak components. 

Two of the main components (Components 1 and 2) 
are identified with the two components detected in the 
much shorter line of sight toward a CMa (Sirius) which is 
not far from e CMa in the sky. One of them (Component 1) 
is the Local Interstellar Cloud (LIC) in which the Sun is 
embedded. For the four elements for which reliable mea- 
surements exist for the two lines of sight (i.e. performed 
with unsaturated lines) we find a constant ratio of 1.5 be- 
tween the column densities of the LIC toward e CMa and 
the column densities of the LIC toward a CMa. 

We derive the neutral hydrogen column density from 
our measurement of O I which is a good tracer of H I. 
Depending on the oxygen abundance we adopt, we find 
a neutral gas column density for the whole line of sight 
between 6. Oil. 2 10 17 cm -2 if we consider that the local 
gas is not depleted (adopting the B stars abundance) and 
9.0±2.0 10 17 cm -2 if we adopt the mean ISM oxygen abun- 
dance of Meyer et al. (1998). With the same two alterna- 
tives, we derive for the LIC a neutral gas column density 
between S.O+o.l 10 17 cirT 2 and 4.4±J;{j 10 17 cm- 2 . 

We estimate the temperatures and electron densities in 
the three main components by combining the information 
of the two ratios N(Cii*)/N(Cii) and N(Mgn)/N(Mgl). 
In particular for the LIC we find n c = 0.12 ± 0.05 cm -3 
and T = 7000±1200 K, both in agreement with previous 
determinations having similar error bars. 

We compare the neutral gas column densities with the 
total (neutral and ionized) gas column densities derived 
from the S 11 measurements, to conclude that Component 3 
is ionized by more than 95 %, that Component 2 is proba- 
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bly neutral but could be as much as 43% ionized and that 
the LIC can be at most 55 % ionized. We conclude from 
these numbers that Component 3 must be located further 
away on the line of sight and is thus almost fully ionized 
by e CMa and that the ionization fraction in the LIC is 
compatible with the gas being ionized by the local EUV 
radiation fields from the hot stars and the cloud interface 
with hot gas. In contrast, it is hard to explain the low state 
of ionization of Component 2 unless it is included within 
the LIC which shields it from the ionizing radiation. 

We detect high ionization species. Sim is detected in 
all clouds but more significantly in the LIC and Compo- 
nent 3, and is probably located in extended layers in the 
outer regions of the clouds. C iv is also detected in the LIC 
and Component 3, but with small velocity offsets from the 
lower ionization species. The derived amount of highly ion- 
ized gas and the derived high temperature are consistent 
with the predictions of Slavin (1989) for a conductive in- 
terface between the Local Cloud and the surrounding hot 
gas from the Local Bubble. 
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